Renormalization group limit cycle for three-stranded DNA 
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We show that there exists an Efimov-like three strand DNA bound state at the duplex melting 
point and it is described by a renormalization group limit cycle. A nonperturbative RG is used to 
obtain this result in a model involving short range pairing only. Our results suggest that Efimov 
physics can be tested in polymeric systems. 
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Consider a three particle quantum system with pair- 
wise short-range potential. Apart from the occurrence of 
the usual three-body bound state, a very special phe- 
nomenon occurs at the critical two-body zero-energy 
state. An infinite number of three-body bound states 
appear though the corresponding potential is not appro- 
priate to bind any two of them; the removal of any one 
of them destroys the bound state. This phenomenon, 
valid for any short range interaction, is known as the 
Efimov effect. The corresponding three-body bound 
states are called Efimov trimers. These states have the 
property that two successive energies, -En's, maintain 
E n+ \/ E n = constant. The size of the trimer is large com- 
pared to the potential range, and so it is a purely quan- 
tum effect [H- Although it was predicted in the context 
of nuclear physics 0, [j], it has now been detected in cold 
atoms Q. 

An ideal DNA consisting of two Gaussian polymers in- 
teracting with native base pairing undergoes a critical 
melting transition where the two strands get detached. 
Maji et.al. recently showed that if, to a double stranded 
DNA at its melting point, a third strand is added, the 
three together would form a bound state instead of re- 
maining critical^]. The existence of a three strand bound 
state has further been verified by real space renormaliza- 
tion group on hierarchical lattices and transfer matrix 
calculations in real space [HQ. That this is an Efimov- 
like effect can be seen by the imaginary time transforma- 
tion of the quantum problem in the path integral formu- 
lation. The paths in quantum mechanics are identified as 
Gaussian polymers and the equal time interaction maps 
onto the native base paring. Such a bound state of a 
triple stranded DNA is called an Efimov DNA. 

In both cases, the special effect is due to a long-range 
attraction generated by critical fluctuations at the tran- 
sition point. For the DNA case, the large fluctuations in 
the bubble sizes at the melting point allow a third strand 
to form bound segments with the other two. The power 
law behavior of the size of a polymer is essential to induce 
a 1 / R 2 interaction between any two chains 0, [B[ . 

Universal aspects of polymers are well understood in 
the renormalization group (RG) approach^. In RG, 
a system is described by length scale dependent run- 
ning parameters, which with increasing length scales are 
expected to reach certain fixed points. Unstable fixed 
points describe transitions while stable fixed points de- 



scribe phases. Both single chain and solutions involving 
many chains are well understood in this framework. The 
purpose of this paper is to show that the triple chain 
bound state at the duplex melting point is of a different 
type. This "few chain problem" is actually described by 
a renormalization group "limit cycle" [HHj. 

Limit cycles are isolated closed trajectories in the pa- 
rameter space to which a running parameter falls into 
when the length scale is changed. These occur if the RG 
flow equation or the /3-function has complex fixed points. 
The appearance of limit cycles invokes log periodicity in 
the corresponding three-body coupling in the polymer 
problem. So they break the continuous scale invariance 
around the two-body fixed point imposing a discrete scal- 
ing symmetry, the hallmark of the Efimov states. 

Another motivation of this paper is to emphasize that a 
three chain polymer model, a three stranded DNA in par- 
ticular, by virtue of mathematical similarities, provides 
an alternative system for Efimov physics. Triplex DNA 
is known to occur in nature. The possibility of recog- 
nizing the bound base pairs of a duplex without opening 
it, by forming Hoogsteen pairs, has the potentiality of 
new types of antibiotics. In addition, H-DNA is a com- 
mon motif formed during many DNA activities where 
there is a stretch of triplex DNA via a strand exchange 
mechanism^. The advantage with our model is that 
the parameters corresponding to the polymers are eas- 
ily tunable by changing temperature, solvent conditions, 
etc. Some important phenomena are getting verified in 
the mathematically analogous low energy condensed mat- 
ter systems; e. g. Majorana fermionsfioj]. the Klein 
paradox [IH. structure formations in early universe [l2|. 
We hope that this work will inspire experimental searches 
of the Efimov effect in polymers. 

The difficulty of the problem can be gauged by consid- 
ering the traditional approach of representing the poly- 
mer Hamiltonian as the Edwards path integrals. In this 
approach, the two chain and the three chain interactions 
have different critical dimensions, namely d = 2 and 
d = 1 respectively. Consequently, in three dimensions, 
one has to tackle irrelevant variables, something outside 
the scope of perturbative RG. 

Our model consists of three directed polymer chains in 
3 + 1 dimensions, where 3 is the dimension of space in 
which monomers live and we assign an extra dimension 
along the contour of a polymer. Two polymers can in- 
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FIG. 1. (Color online) Schematic diagram of strand exchange 
and the equivalent coarse-grained three chain interaction. A 
single strand (solid black line) pairs with the strands (dash- 
dot lines (left)) of a bubble on a duplex (thick line (right)). 
Short green dotted lines indicate base pairings. 



teract only when they are at the same space and length 
coordinate (native base pairing). To avoid difficulties in 
solving the full three chain problem, we consider a simpli- 
fied model that any two of the three chains are in a bound 
state and we do not allow any bending of the bound seg- 
ment. At zero temperature, the two chains form a rigid 
rod-like bound state. At finite temperature, thermal fluc- 
tuation melts locally the bound state to form a number 
of bubbles, or segments of free-chain pair, creating in- 
terfaces. A weight factor 172 is given to each interface 
or Y-fork. The bubbles allow the paired bound state to 
bend. Any unbound chain is taken as free and Gaussian. 
A bubble of two free chains allows a third chain to pair 
with one of them (Fig. [TJ . This is the strand exchange 
mechanism already alluded to. Such an exchange gener- 
ates an effective three chain interaction. In presence of 
large fluctuations in bubbles, we show that the effective 
three chain interaction creates a triplex bound state at 
the critical point of the two-chain melting. 

Our methodology is to find the two chain partition 
function at the melting point and then use this partition 
function to determine the third virial coefficient of the 
three chain problem. Then we calculate the renormal- 
ization group beta function for the effective three-chain 
interaction. Denoting the dimcnsionlcss three-chain cou- 
pling constant as H , our main result is to show that the 
beta function is of the form 



A 9A = ~MH - H )(H - HS), 



(1) 



where A is a real constant and Hq is a complex number. 
Because of this pair of complex fixed points there will be 
a limit cycle behavior, but most importantly, there will 
be a bound state since H — > —00. 
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FIG. 2. Basic building blocks, (a) represents Z(k., s) for a 
Gaussian chain, (b) Zb(k,s) for a two chain bound state, (c) 
a Y-fork representing the interface between a bound pair and 
two open stands. It has a weight 32- 



Let us first define the basic building blocks of the 



model. To avoid the infinite entropy per unit length of 
this continuous model, the unconstrained entropy of a 
single chain is taken as In /x per unit length, and we in- 
troduce an ultraviolet cut-off A. For the two-chain bound 
state, we assign an energy oc e(< 0) per unit length. We 
take fcfiT = 1 where fes is the Boltzmann factor and T 
is the temperature. The free chain and the bound state 
partition functions are given respectively by 



Z(r, N) 
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Zb(iV) = A 3 e 
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We use appropriate powers of A to make other quantities 
dimcnsionlcss. The factor A 3 is used in the definition of 
bound state to make it dimensionally equivalent to the 
free chain. It is convenient to work in the Fourier-Laplace 
(k, s) space where k and s are Fourier and Laplace con- 
jugates to r and N respectively. In (k, s) space the above 
partition functions read 

Z(k,s) = ^ — , Z b (k,s) = T ^ f <f(k).(3) 



s — A 2 log fi 
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Here, k is the Fourier vector conjugate to the end-to- 
end distance. The Gaussian behavior is reflected in the 
average size as measured by the mean squared distance 
(r 2 ) - N. The poles in s of Eq. © for k = at A 2 log/x 
and — eA 2 are the negative free energies of a free chain 
and a bound pair respectively. 
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FIG. 3. (a) The duplex partition function as an infinite 
series of bound pairs and bubbles, (b) Y-fork for a duplex, 
(c) A three chain interaction, (73, involving a free-chain and a 
duplex. 



We are interested in small k behavior. For that we 
treat k as an external variable and drop the 5-function 
for further calculations. The basic building blocks of our 
model are shown in Fig. [2] We allow arbitrary number 
of bubbles. As our partition functions are translationally 
invariant and a forking can take place at any s, we have 
k, s conservation at each vertex point. As arbitrary num- 
ber of bubbles are allowed, the finite temperature bound 
state partition function can be written as an infinite geo- 
metric series as shown in Fig. [3^- We name the black box 
on the l.h.s as a duplex. The duplex partition function, 
obtained by summing the infinite series, is given by, 



Z d (k, S ) = Z b (s) [1 - g 2 2 I Z b (s)Y 
Here the single bubble contribution Iq isfu 
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- '-y/s' + k 2 /A 
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FIG. 4. Diagrammatic representation of the three chain par- 
tition function. This translates into an integral equation in- 
volving interactions to all order. 



in the limit (s' + k 2 /4) — !• where s' = s — 2A 2 log /i. 

We want to concentrate on the events near melting. 
For small s' and k = we find the pole of Z d to be 



At 



Here so' is the free energy difference of the duplex and 
two unbound chains, At is the deviation from the duplex 
melting point and the melting transition is at At = 0. We 
choose g2 as the tuning parameter to get melting. The 
usual scaling for ideal DNA melting can also be recovered 
from Eq. [6j viz., 



so 



C , with £ 



lAtl 



(7) 



Here £ is a diverging length scale for this continuous melt- 
ing. In absence of gi the denaturation transition is first 
order with s = — eA 2 and s = 2A 2 ln/z for the bound and 
the free phases respectively. 

When the system length scale £ diverges for some crit- 
ical value of #2 the system goes to the stable free chain 
phase. The full duplex partition function can now be 
written in terms of £ in the small s' limit as 



Z d (k,s) = (2Tr 2 gl [-r 1 + y/*' + k 2 /4\) 



(8) 



We do not get into details of the two chain melting as 
it is similar to the necklace modelflijj. We note two fea- 
tures. First, there is a diverging length scale £ ~ |At| _1 . 
Second, under a scale change k — > 6 _1 k, the length scale 
and the free energy change as £ — > 6£ and / — > b~ 2 f, re- 
spectively, for any arbitrary b. This is continuous scale 
invariance. 

Now let us consider the three chain case where two 
of them are in duplex state and the third is free. In the 
presence of bubbles we replace Fig. &p by Fig. [3p and call 
it as our two-chain coupling 172 ■ A duplex can dissociate 
into two free chains and if one of them interacts with 
the third free chain they can again form a duplex. As 
a duplex can bend thanks to the presence of bubbles we 
can also have diagram like Fig.[5p. This is our three-chain 
coupling g 3 . 

We denote free chain-duplex interaction vertex by W. 
This is the third virial coefficient which comes from all 
the three chain connected diagrams as shown in Fig. |4] 
Translational invariance suggests that the partition func- 
tion depends only on the duplex-single chain separation 



at the two end points and not their individual positions. 
The total momenta at each end can therefore be taken 
to be zero. 

The purpose of introducing the three-chain interaction 
is to make the virial coefficient independent of the arbi- 
trary cut-off. This can be achieved by doing a momentum 
shell type integration over thin shells which will also give 
the beta function for the 3-chain interaction. 

The evaluation of diagrams of Fig. 0] gives 

W(k,k',s 1 ,s' 1 ,s)=2g 2 Z{k + k',s-s 1 -s' 1 )+g 3 



+2^2 



2m 



At = 2 log/i + e - 4tt 5 2 A 2 . (6) where, 



11 = Z(q, s)Z(k + q, s — Si - s)Z d (-q, s - s) 

xW(q,k',s, s[,s), 

1 2 = Z(q, s)Z d (-q, s - s)W(q, k', s, s[, s) 



(9) 



(10) 
(11) 



We further simplify our model by averaging over 
all angles and assuming the external free chains are 
in relaxed state such that si = A 2 \og fi — k 2 /2, and 
s[ = A 2 log/i — k' 2 /2. The angle averaged partition func- 
tion from Eq. © is then, with s" = s — 3A 2 log/i, 

W(k,k') =g 2 L(s",k,k')+g 3 
+4n^dqq 2 [2g 2 L(s", k, q)+g 3 ] Z d Lq, W(k, q),(12) 

i ti h \ 1 i a + b 2 +c 2 + bc 

and L(a, b,c) = — log -. (13) 

be a + tr + c z — be 

Near melting, s" is very small as there are three chains 
now, and, £ — > oo. So most of the contribution comes 
from the loop diagrams. Neglecting the tree level contri- 
butions and using Eq. ((U in Eq. (fl"2"|). we have 



W{k) 



8k r 



dq 



L(0,k,q) 



H(A) 



A 2 



W(q), (14) 



in terms of the redefined quantities 

W(q) = qW(q), and H(A) = 



.93 

1.9 2 2 



A 2 . 



(15) 



For H = 0, there is no scale in the limit A — > oo, and 
Eq. ([14]) has a solution 



W(k) =Ccos 



*o log £ 



(16) 



where C, A* are constants and so = 1.5036. When H ^ 
and A is finite we can still use this solution as it retains 
its form changing only its constants (l5j. 

We now use a momentum shell technique to get the be- 
havior of H(A). We first integrate over a small shell of ra- 
dius Ae~ dl in Eq. ([T4"|) and then rescale back A — > Ae dl to 
get the following differential equation in the limit l;«A, 



1 

A 



W(q)dq + [1 + H]W{A) = 0. (17) 
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FIG. 5. (Color online) (a) Schematic plot of / vs. T. The 
solid red (dashed black) line is the free energy curve for the 
duplex (triplex) measured from the unbound state. The con- 
tinuous melting is at T = T c . The filled circle is the triplex 
bound state at T c . The triplex melting at T t is first order, 
(b) The corresponding Efimov plot in the K-g plane with line 
types as in (a). 



Using the form of W of Eq. (Q~6|) , we have the solution 
sin (s log 



H(A) 



arctan — 

, so 



sin ( s log 



arctan — 

so 



(18) 



Eqs. (|17p and l|T8j) can be combined to obtain the RG 
flow equation of H as 



A 



dH 
~dA 



l + s 2 



1 + H 2 ) + (1 - s 2 Q )H. (19) 



l+isn 

1—iso 



and Hn 



Eq. (|19p has complex fixed points at Hq 
If we define a new coupling constant £ = ^Zff* > its flow 
gives limit cycle trajectories[l6|. That H is independent 
of A for large values of A and it is log periodic, can easily 
be verified from Eq. (fT5)l . As g 2 ^constant, at the two- 
chain critical point, g% has the same flow equation and 
same form as Eq. (|T9l and Eq. (|T5)) respectively. At the 
points A„ ~ A* (e?) , where n's are integers, 53 runs 
into negative infinity as A is changed [l7j. As 33 can be 
interpreted as three-chain binding energy, wc get Efimov 
states at these points. When A is increased, (73 contin- 
uously decreases to negative infinity and then jumps to 
positive infinity in a log periodic manner. Every jump 
corresponds to a winding around the limit cycle or from 
one Efimov bound state to other. These states are con- 
centrated at the origin and infinite in number. 

So far as the DNA is concerned, the two chain melt- 
ing has a critical behavior with free energy given by 
/ ~ — \T — T c \ 2 for T — > T c — , measured from the un- 
bound state. The spatial length scale £, coming from the 
fluctuations in the bubble sizes, diverges as given by Eq. 
([7]). However the addition of an extra similar strand, de- 
stroys the continuous scale invariance mentioned below 
Eq. ©. Instead of three critical pairs, one gets a fluc- 
tuation induced bound state with a characteristic length 
scale A*. By using the quantum path integral to poly- 
mer mapping, the DNA partition function can be written 
as Z(N) ~ ^2 n exp(—E A 2 l N), where E is the ground 
state energy determined by A*, in units such that EqNA 2 
is dimcnsionlcss. Wc assume that the coefficients do not 



depend too sensitively on n. For N —> 00, the thermody- 
namics is determined by the ground state energy Eq. If 
/ and /at represent the free energy per unit length in the 
long length limit and for finite length N, then one sees 
that 



N 



1 

N 



In 



1 - 



exp(-aN/ aA r) 



exp(-7V/ w ) 



(20) 



valid only for a = cxp(27r/so). 

Fig. [5] shows the equivalence of the free energy curve 
and the conventional Efimov plot in QM. The filled circle 
at T c is the schematic three chain bound state free energy. 
The bound state of size A" 1 , much larger than the hydro- 
gen bond length would melt at a higher temperature at 
T = T t . This free energy curve meets the unbound curve 
at a finite slope indicating a first order transition The 
Efimov DNA is observablcs in the hatched region be- 
tween T c and T t . The equivalent plot for the Efimov case 
is shown in (b). The ground state energy in QM cor- 
responds to the DNA free energy so that in the Efimov 
plot, the y-axis of the wave number (K) in QM becomes 
sign(/) yj/f while the x-axis is the temperature devia- 
tion from the melting point written in terms of the inverse 
duplex length scale, g = sign(T c — T) One recovers 
the square-root behavior of the trimer energy as it ap- 
proaches zero and the duplex curve becoming straight 
line - familiar from the Efimov plot. 

This paper presents an example of the Efimov effect 
which is more amenable to experimental verification. 
Here we show that there exists a triplex state at the crit- 
ical point of duplex melting of DNA. The derivation uses 
nonperturbative RG of a model involving a short range 
pairing only. We start from a zero temperature, or, com- 
pletely bound two-stranded DNA. At nonzero tempera- 
ture, bubbles form in the bound state and a third strand 
can form a duplex with any one of the denatured pair 
or both. The renormalization of the short range pairing 
generates an effective three-chain interaction which is re- 
sponsible for the three-chain bound state. As a result, 
at the critical point of two-chain melting, there exists a 
three-chain bound state, but no two-chain bound state. 
The parameters in the case of polymers are easily tun- 
able, and, therefore, it will be more helpful in detecting 
the Efimov effect experimentally. 
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